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Abstract 



It is proposed that T violation in physics, as well as the masses of electron and 
u, d quarks, arise from a pseudoscalar interaction with a new spin field t(x), 
odd in P and T, but even in C. This interaction contains a factor 275 in the quark 
and lepton Dirac algebra, so that the full Hamiltonian is P, T conserving; but 
by spontaneous symmetry breaking, the new field t(x) has a nonzero expectation 
value < r >^ that breaks P and T symmetry. Oscillations of t(x) about its 
expectation value produce a new particle, the "timeon". The mass of timeon is 
expected to be high because of its flavor-changing properties. 

The main body of the paper is on the low energy phenomenology of the timeon 
model. As we shall show, for the quark system the model gives a compact 3- 
dimensional geometric picture consisting of two elliptic plates and one needle, 
which embodies the ten observables: six quark masses, three Eulerian angles 
#12, #23, #31 and the Jarlskog invariant of the CKM matrix. 

For leptons, we assume that the neutrinos do not have a direct timeon inter- 
action; therefore, the lowest neutrino mass is zero. The timeon interaction with 
charged leptons yields the observed nonzero electron mass, analogous to the up 
and down quark masses. Furthermore, the timeon model for leptons contains two 
fewer theoretical parameters than observables. Thus, there are two testable re- 
lations between the three angles 9 12 , 2 s, 31 and the Jarlskog invariant of the 
neutrino mapping matrix. 

PACS: 12.15. Ff, 11.30.Er 

Key words: timeon, CP and T violation, CKM matrix, neutrino mapping matrix, Jarlskog 
invariant 
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1. Introduction 



We suggest that the observed CP and T violations are due to a new P odd 
and T odd spin zero field r(x), called the timeon field; the same field is 
also responsible for the small masses of u, d quarks, as well as that of the 
electron. Consider first the quark system. Let and be the quark 
states "diagonal" in W ± transitions[l]: 

«d)^ft(T)+w" (i.i) 

and 

fc(T)^ftd)+w + (i.2) 

with i = 1, 2 and 3. The electric charges in units of e are +| for ^(f) and 
— | for These quark states and are, however, not the mass 

eigenstates d, s, b and u, c, £. We assume that the mass Hamiltonians 
for <ft(f) and for are given by 



#T/I = 



01, 02, 03 



T/I 



(G74 + ^7475)1/1 



(1.3) 



with the 3x3 matrices Gyi and F 



'T/I 



92 

v 93 y t/i 

both real and hermitian. The mass 



G Vl - 



;i-4) 



T/I 



matrix Gf/j_ is the same zeroth order mass matrix as Mo(qyi) of Ref. 1, given 
by 

-££77 -a£ a + (3 ) 

in which aj, /?•]•, £f, and q;j_, ^, r/j are all real parameters with ay^ 
and to be positive. It can be readily verified that the determinants 

|<3 T | = Knl = o. (1.5) 

Thus, the lowest eigenvalues of and Gj, are both zero. These two real 

symmetric matrices can be diagonalized by real, orthogonal matrices {U^)q 
and {Ui)o, with 

/ N 



(t/ T )jG T (t/ T ) = 



m (c) 
m (t) y 



(1.6) 
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and 

/ \ 

(^1)^(^)0= mo(s) (1.7) 

V m (b) ) 

where the nonzero eigenvalues are the zeroth order masses of c, t and s, b 
quarks, with 

m (c)=/? T [l + rj\{l + $)], (1.8) 

m (*) = a T (l + £ 2 )+/? T , (1.9) 

mo(s)=A[l + 77 2 (l + ^ 2 )] (1.10) 

and 

m (6) = « i (l + ^ 2 ) + /3 i . (1.11) 

Thus, Gj and Gj. can be each represented by an ellipse of minor and major 
axes given by rao(c) and mo(t) for f and likewise mo(s) and mo (6) for j. 

The orientations of these two elliptic plates are determined by their eigen- 
states. As in Ref. 1, we define four real angular variables 0j ar| d @h H 
by 

^ = tan i? £ T = tan0 T 
r/l = tan 6*| cos 0| and r/j = tan 6*| cos 0f. (1-12) 
The eigenstates of Gf are 

/ cos6*t \ 



sin cos 0| 
\^ sin #| sin 0j y 



with eigenvalue 0, 



( 



— sin #f 

COS Q\ COS 0| 

cos B\ sin 0| 



with eigenvalue mo(c) 



(1.13) 



(1.14) 



and 





- sin 0| 

COS 0| 



with eigenvalue mo(t). 



(1.15) 
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Correspondingly, the eigenstates of Gj are 

( cos 6 1 \ 
— sin 9 1 cos 4>i 



\ - sin 9 i sin 



with eigenvalue 0, 



(1.16) 



PI = 



( sin^ \ 

cos 6 1 cos 4>i 
cos 9 l sin 0| J 



with eigenvalue rao(s), 



(1.17) 



and 



Pi = 







V 



■ sin (p[ 

COS 01 



with eigenvalue mo(b). 



(1.18) 



We note that by changing 0f, 0| to — 4>i the unit vectors e|, p|, Pf of 
(1.13)-(1.15) become ej_, pj_ and Pj_ of (1.16)-(1.18). Here the signs of 6^ and 
are chosen so that the sign convention of the particle data group's CKM 
matrix agrees with both 9^ and 9^ being positive, as we shall see. In terms of 
these eigenstates, the 3x3 real unitary matrices (£/f)o and (£/j,)o of (1.6)-(1.7) 
are given by 

(t/ T ) = (e T , Ph P T ). (1.19) 

and 

(Ui)o = (e b P b Pi) (1.20) 

Thus, in the absence of the iFy^j^ term in (1.3), the corresponding CKM 
matrix in this approximation is given by 



(Uckm)o = p/ t )J(^)o = 



cos 9 1 cos 6^ 

— sin 9 1 sin 0f cos 

— cos 6*i sin 

— sin 6*i cos cos 

sin #i sin 



sin 1 cos 6*i sin 0-f sin 4> 

+ cos 6*i sin 6*i cos0 

— sin 9 1 sin <9j cos 0f sin 
+ cos 9 1 cos 6*i cos0 



\ 



— cos 9 1 sin (f) 



COS0 y 



(1.21) 
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in which 

We assume that T violation and the small masses of u, d quarks are due 
to the new 

zF 74 75 (1-23) 

term in (1.3), with 

F ] = F i = F = r q ff (1.24) 

in which r q is a real constant and / a 3 dimensional unit vector represented by 
its 3 x 1 real column matrix. Graphically, we can visualize Gf and G[ as two 
elliptic plates mentioned above, and Fy^ as a single needle of length r q and 
direction /, as shown in Figure 1. 

We note that the time-reversal operation T in quantum mechanics involves 
a complex conjugation operation changing the factor % to —i. Since the entirety 
of classical mechanics can be formulated with real numbers, the presence of i 
in quantum mechanics is necessitated by the commutation or anticommutation 
relation between operators, such as that between 74 and 75. This led us to 
postulate (1.23) as the source of T violation. The specific form given by 
(1.23)-(1.24) can be due to the spontaneous symmetry breaking of a new T 
odd, P odd and CP odd, spin field r(x), which has a vacuum expectation 
value given by 

< T{X) > vac = T q ^ 0. (1.25) 

While the general characteristics of spontaneous time reversal symmetry break- 
ing models have been discussed in the literature[2], one of the new features of 
the present model is to connect such symmetry breaking with the smallness of 
the light quark and electron masses. 

In Section 2, we begin with a general 3 x 3 T, P and CP odd mass matrix 
of the form 

£74 + ^7475 (1-26) 

with Q and T both real and hermitian, then derive some useful properties of its 
eigenvalues and eigenvectors. In Sections 3 and 4, we summarize the analysis 
of how in (1.24), the length r q and the direction / of the needle are related to 
the light quark masses and the Jarlskog invariant[3] J of the CKM matrix[4,5]. 
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As we shall see, this leads to 



r q ^ 33MeV, (1.27) 

m u S r g (/e T ) 2 (1.28) 

and 

^ = r g (/ ei ) 2 (1.29) 

with e-f and ej_ given by (1.13) and (1.16). 

An interesting feature of the model is: the implicit assumption that the 
constant r q might be due to the spontaneous symmetry breaking of a new 
type of T odd and CP odd, spin field r{x), which has a vacuum expectation 
value given by (1.25). As an example, we may assume that the Lagrangian 
density of t(x) is given by 

1 / dr \ 2 

~V(t) (1-30) 

with 



V(r) = ~\^\t\ - \t 2 ) (1.31) 

in which the (renormalized) value of A is positive. This then yields (1.25). 
Expanding V(r) around its equilibrium value r = r q , we have 

V(t) = - \t\ + \ml(T - T q f + 0\(t - r,) 3 ] (1.32) 

with 

m T = (2A)^r g , (1.33) 

the mass of this new T violating, C violating and CP violating quantum, called 
timeon. 

The interaction between r(x) and the quark field might be obtained by 
replacing the F = r q ff factor of (1.24) with 

F = r(x)ff. (1.34) 

Because of the flavor-changing property of the timeon field [6] , its mass (if 
it exists) must be quite high. A full analysis of this interesting possibility 
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lies outside the scope of this paper. Here, we concentrate on the low energy 
phenomenology of the timeon model. 

In the application to quarks, there are ten measurable parameters in Hy^ 
given by (1.3). These consist of 3 angles 

9 h 6 i and = T -^ (1.35) 

of (1.21) and 4 zeroth order masses 

mo(c), rao(t), mo(s) and mo (6) (1.36) 

given by (1.8)-(1.11) in the description of and In addition, the timeon 
term (1.24) contains 3 parameters: 

T q and two angles in /. (1-37) 

(The angle 0| + 0j_ is an unphysical gauge parameter.) These ten theoretical 
parameters account for ten observables: six quark masses, three Eulerian angles 

#12, #23, #31 (1.38) 

and the T-violating phase factor 

e lS (1.39) 

in the CKM matrix. The timeon model provides a simple picture, combining 
these ten observables into a single compact geometric structure. 

In Section 5, we extend the timeon model to leptons. Similar to quarks, 
there are also 10 observables: six leptonic masses and four angles (as in (1.38) 
and (1.39)) of the neutrino mapping matrix. However, unlike the ] and J, 
quarks, the mass scales of the neutrinos are much smaller than those of the 
charged leptons. Thus, it seems reasonable to explore the interesting possibility 
that the timeon term is absent in the neutrino sector. In this sense, the neutrino 
sector may be regarded as more " primeval" , with its lowest neutrino mass zero. 
As we shall discuss, the corresponding timeon model for leptons gives electron 
a mass and the number of parameters in the theory can be reduced to only 8. 
Thus, there are 2 testable relations between the 10 observables. 
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2. Eigenvalues and Eigenvectors 
2.1 General Formulation 



We begin with (1.26) and write the corresponding mass Hamiltonian as 

H = ^(^74 + ^7475)^ (2-1) 

where Q and T are both 3x3 real hermitian matrices. Resolve the Dirac field 
operator ip into a sum of left-handed and right-handed components: 

^P = L + R (2.2) 

with 

L = i(l + 75 )^ and fl = i(l- 75 W>. (2.3) 
Thus, (2.1) becomes 

H = L\Q - iT) lA R + R\G + iT) lA L. (2.4) 

Define 

M = Q - iT, (2.5) 

M^ = Q + iT (2.6) 

and their product 

M 2 = MM^ = {Q- %T){Q + iT). (2.7) 

Since A4 2 is hermitian, it can be diagonalized by a unitary transformation. We 
write 

V ] L MM^V L = m 2 D = diagonal (2.8) 

with 

vlv L = 1. 

Multiply (2.8) on the right by m^ 1 , we see that by defining 

V R = M^V L m D \ (2.9) 

we have 

VlMV R = m D ; (2.10) 



furthermore, Vr is also the unitary matrix that can diagonalize the correspond- 
ing M f M; i.e., 

V ] R M^MV R = m 2 D . (2.11) 

(Note that M 2 = MM\ but M can be quite different from M or .) 

For application, the matrices Q and ^ can be either Gp Ff or Fj_ of 
(1.3). As in (1.13)-(1.15), we write the eigenstates of Q as column matrices 



e, p and P. 

Likewise, write the unit column matrix / in (1.24) as 



(2.12) 



/ = 



cos a 

sin a cos b 
sin a sin b 



Thus, 



Q = uee + fipp + mPP 
and correspondingly we may set 

T = Tff, 



(2.13) 



(2.14) 



(2.15) 



with is, fi, m, r all real constants. When r = r q and v = 0, T becomes F of 
(1.24) and Q can be either G-\ or Gj_ of (1.4). For the moment, we retain the 
eigenvalue v in (2.14) for the formal symmetry of some of the mathematical 
expressions (as in (2.23) below), even though v = when we discuss physical 
applications of our model. 

As in (1.19) and (1.20), we define a real unitary matrix Uq whose columns 
are e, p and P of (2.12); i.e., 



U Q = (e p P). 
The matrix Uq diagonalizes Q, with 



(2.16) 



G' = UoGUo = 



I 



v \ 

fJL 

m j 



(2.17) 
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It also transforms T into 



T' = U TU = rf'f 



where 



(fe \ 

fp 
\fp ) 



with 

As before, 



fe = e/, f p = pf, fp = Pf. 



2.18) 
2.19) 
2.20) 

/e+/p+/p = l- ( 2 - 21 ) 

By using (2.17)-(2.19), we find that the same Uq also transforms the matrix 
M 2 into 

{M'f = U M 2 U = {G'f + {T'f + i[Q' 
which is given by 

/ v 2 + r 2 f 2 t[t - v)]f e fp T\r-i{m-v 

T[T +i(fJL - V)]f e fp l^ + T 2 / 2 T[T-i(m-fJL 

t[t + i(m - v)]fpf e t[t + i(m - ti)]f p f P m 2 + r 2 fj> 



{M'f = 



For our applications, we are only interested in the case v = 0. Define 

M = \im{M'f 

and let A 2 , A 2 ,, A 2 , be the eigenvalues of N. From (2.23) and (2.24) we 

A 2 + A 2 + A3 = m 2 + ^i 2 + r 2 

and 

A?A|A§ = PV| = SfU?m\ 
These eigenvalues are also the solution A 2 of the cubic equation 



\M - A 2 | = \H\ + AX 2 + BX A - A 6 = 



where 



^4 = V™ 2 - r 2 [m 2 (l - / 2 ) 2 + ^ 2 (1 - f 2 f + 2m^f 2 f 2 ] 



2.22) 

]f P fe \ 
}fpfp 



2.23) 

2.24) 

have 

2.25) 

2.26) 
2.27) 

2.28) 
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and 

B = m 2 + fi 2 + r 2 . (2.29) 

In the limit r — > 0, we see from (2.25)-(2.29) that the two heavier masses be- 
come fi and m, while the lightest mass is proportional to r. (Readers who are 
only interested in perturbative solutions are encouraged to move on to Section 
3 directly.) 



2.2 Some Useful Expressions 



It is convenient to arrange the three eigenvalues X 2 , X 2 ,, A3 in an ascending 
order, each with a new subscript: 

A 2 < A 2 < A 2 , (2.30) 

Write 

Ei = X 2 (2.31) 

with 

i = s, I and L. (2.32) 

(Here, the letters are s for small, / for large and L for very large.) Let ipi be 
the corresponding eigenstate defined by 

■Afipi = E^i. (2.33) 
Using (2.23)-(2.24), we can express J\f in the form 

•W = ( \ "J ) (2.34) 



where 



and 



+ 2 



U _ ( r2 fe T(T ~ ifl)f e fp \ (9 

lT(T + We/ P ^+r% 2 J' ( 35 j 

n 2 = m 2 + r 2 /|>, (2.37) 
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Correspondingly, each eigenvector ipi of J\f can be written as 

4„ =(J) (2.38) 

with ^ a 2 x 1 column matrix and q a constant. From (2.33), (2.34) and 
(2.38), we have 

hfo + (na)x = Ei<j>i (2.39) 

and 

n X ] & + n 2 Cl = Em. (2.40) 

From (2.39), it follows that 

0, = {Ei - h)- l nc lX i (2.41) 

and likewise from (2.40), 

d = ^ - n 2 )- W^- (2-42) 
Substituting (2.41) to (2.40) we find 

m? = £?i = n 2 [l + X f K 2 - /i)-^]- (2.43) 
Likewise, (2.39) and (2.42) lead to 

2 

(ft " ^—ErXX^i = E i( f>i. (2.44) 
n l — hi 

Both (2.43) and (2.44) are valid for all three solutions i = s, I and L. From 
(2.44), we see that once the eigenvalue E^ is known, the determination of the 
corresponding three-dimensional eigenvector ipi reduces to the much simpler 
two-dimensional spinor equation (2.44). This fact will be of use in the later 
sections on leptons. 

Remarks It is well known that a mass matrix of the form (2.4) can also be 
written as a single term 

^ f M74^, (2.45) 

but with the hermitian matrix M complex. However, the reality conditions 
of Q and T lead to a specific form of M, and that specific form may appear 
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"unnatural" without the "timeon" picture, Further discussions will be given in 
Appendix A. 



2.3 Determination of Eigenvalues 



For quarks, the top mass m t is much larger than m c and m u , so is the 
bottom mass m& >> m s and m^. Likewise for charged leptons, m T is >> m M 
and m e . Thus, in the notations of (2.32), for % = L the largest mass m L 
satisfies 

ttil >> mi > m s . (2.46) 

Eq.(2.43) gives a convenient route to express mi in terms of the parameters 
in J\f, given by (2.34)-(2.37). In this case, we can regard the parameter m in 
(2.23) as satisfying 

m = 0(m L )»id and r. (2.47) 
Define the parameter e through 



2 2 , 2 i 2 .f 2 i 

m L = n +e = m +r/ P + e. 
Eq.(2.43) in the case z = L gives 



2 2 t 

m x = n + x 



n 2 + e — h 



X 



which leads to 



e = X 



t. 



n 2 + e — h 



= xUl- 



x = x 



h - e\-i 



h — e 



= x ] x + x ] — ^x + x f 



77/ 



x + 



Thus, we have the expansion 

i 2 L = m 2 + r 2 f 2 P + X f X + - Wto - (X f x) 2 ] + 0(m- 4 ): 



i.e., on account of (2.35)-(2.37) 



mi = m l + T 2 f P l + ft +(1 



iiV 



m' p 



+ 



2 2 
T jl 

m 2 



;2.48) 



(2.49) 



(2.50) 



(2.51) 
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To derive the corresponding expansions for the smaller masses m 2 and mf, 
define 

S = m 2 + mf 

and (2.52) 

V — m 2 mf. 

Thus, 

m 2 = \[S + {S 2 - 4P)5] 
and (2.53) 

m 2 = ±[<S-(<S 2 -47^]. 

From (2.25)-(2.26), we have 

5 = m 2 + ^i 2 + r 2 - m\ 
and (2.54) 

P = T 2 f^ 2 m 2 /m 2 L . 

Combining these with (2.51), we have the expressions for mf and m 2 . In the 
limit m — > oo, but r comparable to //, (2.54) becomes 

<S^ 2 + r 2 (/ 2 + / 2 ) 2 
and (2.55) 

V - r 2 / e V- 

If in addition r — > 0, then we have 



m 



m 2 , m 2 — > ^i 2 and m 2 — ► r 2 / e 4 . (2.56) 
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3. Perturbative Solution and Jarlskog Invariant 



3.1 Perturbation Series 



In this and the following sections we return to the mass Hamiltonian (1.3) 
for quarks and calculate its eigenstates by using G74 as the zeroth order Hamil- 
tonian and 2^7475 as the perturbation. From the discussions given in the last 
section, we see that this is identical to the problem of finding the eigenstates 
of N regarding r as the small parameter. Using (2.23)-(2.24), we may write 



N = 



( 
11 2 








m z 



J 



+ M + 0(r 2 ) 



(3.1) 



with 



M = 





¥fefp 



V 



(3.2) 



-v/e/j -imfpf e N 

-i(m - n)f p fp 

imfpfe i(m-fj)f p fp 

To first order in M\, the eigenstates of J\f can be readily obtained. For ap- 
plications to physical quarks, we need only to identify that f e , f p and fp are 
replaced by 

(/e)t/| = feyi, 

(fph/i = fPyi ( 3 - 3 ) 

and 

(fp)vi = f p yi- 

Likewise, neglecting 0(r 2 ) corrections, we can relate fj 2 and m 2 to (quark mass) 2 

by 



2 2 

^ = m c , 
2 2 



m 2 b 



and set 



(3.4) 
(3.5) 
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Thus, to 0(r q ), in the f sector the state vectors of u, c, £ are related to those 

of 6|, p-f and P| by 



with 



and 



c 



(€ t |«) (p T |tt) (p t N^ 
(e T |c) (p T |c) (P T |c) 

(e T |f) CPtl*) ( P tl*) 



(e T | W ) = l + 0(r a 2 ) : 



(et|c) = -*-^(/e/p)T + 0(r^) : 

/ / br- 



(^\t) = -i±(fjph+0(r^, 

Tilt 



(PTW = -i-r(/e/p)T + 0(r, 2 ), 



m 



b T |c) = l + 0(r 2 ) : 



(Ptl*) = 



m ( + m, 



■(/p/p)t + O(t'), 



(F,\u) = -i-HfJph + 0(Tj), 

Tilt 



nit + rn< 



(fpfph + 0(T 2 q ) 



(P T |t) = l + 0(r 4 2 ). 



Likewise, for the [ sector, we may write 



(d\ 


( 


s 






\ 



>ll s ) (Pi\s) (Pi\s 



\ 


<*l ) 




Pi 


J 





with 



(e l \d) = l + 0(r*) : 



^i\s) = -i±(fef P )i + 0(ri) t 

Tllg 



(3.6) 

(3.7) 
(3.8) 

(3.9) 

(3.10) 

(3.11) 
(3.12) 

(3.13) 
(3.14) 

(3.15) 



(3.16) 

(3.17) 
(3.18) 
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(zi\b) = -i±(fefp)l + 0(r), 
nib 



T, 



(Pl\d) = -i^(fef P )i + 0(r) : 



m 



(pj 5 ) = l + 0(r_ 2 ), 



(Pi\b) = 



rrib + m. 



{f P fp)i + 0(r"), 



(P i \d) = -i-^(fJ P ) l + 0(^), 

Till) 



(Pl\s) = 



— %- 



TTlb + ITl; 



ifpfp)l + Otf) 



and 



(P i |6) = l + 0(r, 2 ). 



3.2 Jarlskog Invariant 
Write the CKM matrix as 

UcKM = 



U21 U22 U23 
U31 U32 U33 



Following Jarlskog[3], we introduce 

Si = U^Un, S 2 = U21U22, S3 = Ul x Uz2 

and define 

J = ImSlS 2 . 

By using (3.27) we see that 

J = Im 

Because of unitarity of the CKM matrix, 

Si + S 2 + S3 = 0. 



(UuU22)(Ul 2 W 2l ) 
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Therefore, J is equal to twice the area of the triangle whose sides are Si, S 2 
and S3. Furthermore, from the explicit form of J given by (3.29), we see that 
J is symmetric with respect to the interchange between the row and column 
indices of the CKM matrix. It follows then in deriving J, we may use the 
elements of any two columns and of any two rows of the CKM matrix. 
It is convenient to denote {Uckm)o of (1.21) simply as V, with 

( Vn Vu V 13 \ 



V — (Uckm)o = 



V21 V22 V23 
V31 V32 V33 



(3.31) 



In terms of the state vectors e^, pp Pf and ej_, Pj_ of (1.13)-(1.18), we 
can also write V as 

( ( e tl e l) ( e tbl) ( e tl p j) \ 

^ = (Ptki) (PtlPi) (Ptln) • (3-32) 

u p ti e i) ( p th) (^in)y 

Likewise, the CKM matrix is given by 

/ (u\d) (u\s) (u\b) ^ 
Uckm = (c\d) (c|s) (c|6) 
^ (f|s) (t|6) J 

= V + ir q W + 0(r 2 q ) (3.33) 

where the matrix elements of W are derived from (3.7)-(3.15) and (3.17)- 
(3.25). Using the perturbative solution of Sec. 3.1, we can readily express 
the matrix elements of Uckm in terms of the corresponding ones of V. The 
Jarlskog invariant can then be evaluated by using (3.29). 

As will be shown in Appendix B, the result, accurate to the first power of 



T q , IS 



J = T Q 



ifJph A + (/e/p)| ^ + 



m s 

(fefph 



rrir 



m s + mfe 
(/p/p)t 



Pi 



m c + mt 



-P 



where 



A s = -^13^23^33 = 



-2 • 10 



(3.34) 



(3.35) 
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A b = -V 12 V 22 V 32 8.8 • 10- 3 , (3.36) 

B i = -V n V 2l V 31 = 1.10 • KT 3 , (3.37) 

A c = V 31 V 32 V 33 -2 • 10- 4 , (3.38) 

A t = V 21 V 22 V 23 = -8.8 • 1(T 3 (3.39) 

and 

£ T = V U V 12 V 13 ^ 1.10 • 10~ 3 . (3.40) 

From the definition (3.29) and (3.34), these coefficients A s , • • • , are 
all products of four factors of Vij. As will be shown in Appendix B, because V 
is a real orthogonal matrix, these quartic products can all be reduced to triple 
products given by (3.35)-(3.40). Since m c » m s and m t » rrib, we can, as 
an approximation, neglect the terms related to the up sector in (3.34). 



4. Determination of r q and / 
4.1 A Special coordinate system 

For the f quarks, the parameters Ai, A2 and A3 in (2.25)-(2.29) are related 
to the quark masses by 

Ai = m u , X 2 = m c and A3 = m t . (4.1) 

Likewise, f £ is 

(/«)t = hi (4-2) 
with e| given by (1.13) and / the unit directional vector of (1.24). We work 
to the lowest order in r q . From (2.26), setting A1A2A3 = m u jim we have 

m u = r,(/e T ) 2 . (4.3) 

Likewise, for the [ quarks 

m d = r 9 (/ ei ) 2 . (4.4) 
It is convenient to introduce a special coordinate system in which 



(4.5) 





(1\ 




/ 


cos6* c ^ 


H = 





and 6| = 




— sin 6 C 








V 


) 
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Since pi and Pj_ are both _L ej_, we may write 





f ^ 




r \ 


Pi = 


— cos 7 


and P| = 


— sin 7 




v sin 7 j 




^ -COS7 ) 



(Uckm)o = 



(4.6) 



Furthermore, we shall set the zeroth order CKM matrix (Uckm)o of (1.21) to 
be 

/ (etki) (e T bi) (e T ln) \ 

W e i) Otbi) btl p i) 
v (Ak) (Ah) (AIA) y 

' .974 .227 5 • 10~ 3 \ 
.227 .973 .04 
5-10- 3 -.04 .999 

Thus, with the same accuracy of O(10 3 ), the Cabibbo angle 6 C is given by 



V 



+ O(l-10~ 3 ) 



(4.7) 



cos# r = .974 and sin# r = .227. 



(4.8) 



Likewise, from (e-f |Pj_) = 5-10 3 in (4.7), in accordance with (4.5), (4.6) and 

(e T |P i )=sin^sin 7 , (4.9) 

we find 

sin 7 = 2.2-Kr 2 , (4.10) 

which together with (4.5) and (4.6) give the coordinate system defined by 
(ej_, pj_,Pj_). Eq.(4.7) then, in turn, determines the corresponding coordinate 
system (e-f , p-f, Pf). 

Next, we shall determine the parameters r q and the directional angles a and 
(3 of the unit vector 

( sin a cos (3 ^ 
sin a sin j3 
cos a 



f = 



\ 



(4.11) 



in the coordinate system defined by (4.5)-(4.6). 
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4.2 Determination of (3 

From (4.5) and (4.11), we have 

/ej_ = sin a cos (3 (4. 12) 

and 

je| = sin a cos((3 + C ). (4-13) 
Thus, on account of (4.3) and (4.4), 

cos 2 (/5 + (9 c ) m u 
cos 2 (3 rrid 

and therefore 



(4.14) 



cos (3 ^irid 

assuming 



cos(/? + # c ) = ± (™tA*_ (415) 



(4.16) 



m d 2' 
we find two solutions for (3: 

/3 = 48°50' (4.17) 



or 

/?^82°20'. (4.18) 
4.3 Determination of a and r q 

We shall first determine the parameter a by using the Jarlskog invariant 

J = 3.08 • 1CT 5 . (4.19) 

Define 

F = 10 2 Jm b /r q . (4.20) 
From (4.4), (4.5) and (4.11), we have 

Trid = T q sin 2 a cos 2 (3 (4-21) 

and therefore 

F = 10 2 J(m b /m d ) sin 2 a cos 2 (4.22) 
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For definiteness, we shall set the various quark masses as 

m d ^ 5MeV, m u ^ 2.5MeV 

m s = 9bMeV, m c = 1.25GeV 
m b ^ 4.2GeV and m t ^ 17bGeV, (4.23) 
consistent with the Particle Data Group values[7]. Thus, (4.22) becomes 

F(a, (3) = 2.6 sin 2 a cos 2 (3, (4.24) 

On the other hand, from (3.34) and by using the numerical values for 
A s , Ab, • • • , #t of (3.35)-(3.40) together with the various quark masses 
given above, the same function F(a, (3) is also 

F(a, /?) * -0.88(/ e /J i - 0.067 (/ e /„) T 

+ 0.88(7^ - 0.021 (/ e / P ) T 

+ 0.H(/ p /p)| +0.0026(/ p / P ) T . (4.25) 

As an approximation, we may neglect the contributions of the f sector. Com- 
bining (4.24) with (4.25), we find 

2.6 sin 2 « cos 2 /? = 0.88/ e J/ A - / P J + 0Al(f p f P ) l (4.26) 

with 

A, = hi (4-27) 

given by (4.12), 

fpi = fPl = — sin a sin (3 cos 7 + cos a sin 7 (4.28) 

and 

fp l = /P| = — sin a sin (3 sin 7 — cos a cos 7. (4.29) 
By using 7 from (4.10), 

/? ^ 48° 50' 

from (4.17), we find 

a ^ -36° 10'. (4.30) 
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From (4.4) and (4.12), we have 

Tq = fTld/ ( s i n2 OL cos2 P) ■ 



(4-31) 



The above values of a, f3 and rrid = bMeV give 

T q = SSMeV. (4.32) 

On the other hand, the alternative solution (3 = 82° 20' of (4.18) leads to 
a much larger value r q ~ 5.5GeV. Such a large value invalidates the small 
approximation. Thus, we focus only on the solution (4.32) in this paper. 



5. Applications to Leptons 

The mapping matrix for leptons presents us with a different quantitative picture 
from the CKM matrix for quarks. 

In the CKM matrix, the only off-diagonal elements of a relatively significant 
size are those associated with the Cabibbo angle 6 C , which is also not large. 
In the lepton mapping matrix Ui, except for its T-violating element ((£//)i3 
in the standard form), all matrix elements are not small. Furthermore, within 
the present ~ la accuracy, Ui is given by the Harrison-Perkins-Scott (HPS) 
form [8] 

This leads us to propose the following simple Timeon model for leptons. 



(Ui)o = 



(5.1) 



5.1 Mass Matrices for Leptons 

As in (1.1)-(1.2), we define ^(|) and ^(J,) to be the lepton states "diagonal" 
in W ± transitions, so that 

«t) Mi) + w+ 

and (5.2) 
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with % = 1, 2 and 3. Their electric charges in units of e are for ^(|) and — 1 
for ^(|). However, these ^(|) and ^(j) are not the mass eigenstates v\, z/ 2 , ^3 
and e, ^i, 



r. Their mass Hamiltonians and are given by 



H 



T/l 



= (/{, it 4) m (C(/)74 + ^(/)7475) T /l 



V /a / T/l 



(5.3) 



We choose and Gj so that if F^n were to vanish, the lepton mapping matrix 
would be (5.1) exactly and, in addition, the lowest lepton mass in either f or 
I would be zero. We further simplify our model by assuming that the timeon 
field does not couple to neutrinos; i.e., 



F(/) T = 0. 



(5.4) 



Thus, the lightest neutrino is massless. All departures of the lepton mapping 
matrix from the HPS form (5.1) stem from the single timeon term 



F(l)i = TiVV 



(5.5) 



for the charged leptons, with v a real 3 dimensional unit vector like / of (1.24) 
for quarks. The same timeon term (5.5) also gives the electron mass. Thus, all 
departures from (5.1) would come from the charged leptons. A consequence 
of the assumption F(l)^ = is that the predictions of the model will now 
depend on the choice (formerly arbitrary) of the "hidden" bases I2, h)yi 
of (5.2). As illustrated in Figure 2, we have chosen this basis to recognize the 
factorization of (5.1) with 

(Ui)o = W)JW)o 



in which 



/ 1 






r 

2 



0\ 



Vf y/l) 



(5.6) 



(5.7) 



and 



Wo = 



' vf o\ 
001) 



(5.8) 
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It appears natural to associate the left-hand factor in the product (5.6) with 
the charged leptons and the right-hand factor with the neutrinos; i.e., through 
(Vj)o of (5.7), we identify 

(h)i = e (5.9) 

and through (V T ) of (5.8), 

(«T = v*. (5-10) 

In (5.9), eo denotes the zeroth order electron state (i.e., without its timeon 
correction). It is convenient to consider (5.6) as a product of two consecutive 
rigid body rotations with eo and v% as their respective fixed axes of rotations 
of the same rigid body. In this picture, we may identify 



(5.11) 



A comparison between (5.9) and (5.10) suggests that like eo, ^3 is also massless. 
(There is an alternative possibility with v\ massless, as will also be discussed 
below.) 

5.2 Analysis of Mass Matrices 

We begin with of (5.3) for the charged leptons by assuming 

/ \ 



G(l)l = 



aj_ + 6j — aj_ 
V ~ a i a l + b l ) 



(5.12) 



with aj_, 6j_ both positive. Using (5.7) we can diagonalize G(l)i through 

( \ 

W)oC(OiW)o = 



(5.13) 



/ 



64. 
\ 2a l + b l 

and leads to the zeroth order mass of e to be 0, those of mu and tau to be 

V = bl 



and 



(5.14) 



m = 2ai + 6j_. 
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The physical masses of e, \i and r can then be obtained by using (2.23)-(2.29), 
with the parameters r and / replaced by 77 and v of the leptonic timeon term 
(5.5). 

In the t sector, in accordance with (5.4) the neutrino masses are only due 
to the term, with the lowest neutrino mass zero. There are two such 

possibilities, of which the first is 

(i) 



G(0 T = 



|a-f- a-f + 6j 



V 



(5.15) 



/ 





where a-f and 6j are both positive. The masses of v\, v<i and are given by 



mi = 6|, m 2 = -a T + 6 T 



and 



m3 = 0. 

In this case, we would have from the present experimental data[9,10] 

\{m\ + ml) ^ 2.39 x lO^eF 2 

and 

m\-m\^ 7.67 x lO^eF 2 . 
(ii) The second possibility is, 

( i a t \/i a T ^ 

G(Z) T = 



2 a T 




a-f 
6 T J 



which leads to 
and 



m\ = 0, 



In this case, instead of (5.18), we have 
and 



m\ = 7.67 x lO- 5 eV 2 



m\ ^ 2.43 x 10- 3 eV 2 . 
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(5.16) 
(5.17) 

(5.18) 



(5.19) 



(5.20) 



(5.21) 



In either case, (i) or (ii) the mass matrix (5.15) or (5.19) can be diagonalized 
with the same real unitary matrix (Vj-)o of (5.8) so that 

( m x \ 



(y T )jG(/) T (F T )o = 



(5.22) 



ra 2 
V ra 3 j 

and the zeroth order neutrino mapping matrix (5.6) is the HPS form (5.1). 

In Fig. 2, we give a graphic illustration of these two rotations (Vj,)o and 
(V|)o- In either (i) or (ii), the parameters in the neutrino sector are determined. 
The remaining parameters are 

/i, m (5.23) 

the two mass parameters of (5.14) for the charged leptons, and the three 
parameters 

T\ and two angular variables (5.24) 

that characterize the unit vector v of the timeon term (5.5). On the other hand, 
these five parameters in (5.23)-(5.24) should account for seven observables: 
the three charged lepton masses 



and the four parameters 
and 



ra e , m M , m T (5.25) 

012, #23, 031 (5.26) 



the Jarlskog invariant Ji (5.27) 

of the lepton mapping matrix. Hence, the model defined in this section pre- 
dicts two relations between these seven observables (5.25)-(5.27), as we shall 
discuss. At present, the existing knowledge of the neutrino mapping matrix 
agrees with the HPS form (5.1) to « 10%. Anticipating that future experi- 
ments may improve the accuracy to 1% level, we shall calculate the elements 
of the lepton mapping matrix to all orders in ri/fi, first order in ri/m, but 
neglecting all 1/m 2 corrections. 

5.3 Statevectors e, ji and r 
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Consider the ^Q) sector. From (5.13)-(5.14), we have 

/ \ 



W)^(OiW)o = 



ii 
\0 m ) 



(5.28) 



As in (2.12) and (2.14), we can express G(l)i in terms of its eigenvalues 0, fi, 
m and their corresponding eigenvectors e, p, P as 



G(l) i = fipp + mPP. 



(5.29) 



Under the same transformation (Vj)o, the T odd timeon term F{l)i of (5.5) 
becomes 

(n)^(Oi(n)o = T~if'f' 

with /' related to the unit vector t> of (5.5) by 



(5.30) 



As in (2.19), write 



( fe \ 

fp 

v fp J 



(5.31) 



(5.32) 



where in place of (2.20) we have 

fe = ev, f p = pv, fp = Pv. 



(5.33) 



Equate the matrices Q and T of (2.1) with G(l)i and the matrix A/" of 

(2.24) becomes 

N={V i )lM\V i ) Q (5.34) 
with M 2 of (2.7) given now by 



M 2 = [G(l) i -iF(l) i ][G(l) i + iF(l) i ]. 



(5.35) 



For applications to charged leptons, we identify the subscripts s, / and L of 
(2.32) in Section 2.2 with e, fi and r. Thus, (2.31) and (2.32) become 



Es = XI = mi 



(5.36) 
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and 

i = e, ji and r. (5.37) 

These lepton masses and their corresponding eigenvectors ip e , ip^ and ip T can 
be readily obtained by using results derived in Sections 2.2 and 2.3, as we shall 
see. 

For i = e and fi, the eigenfunction ipi can be written as, in accordance with 
(2.38), 

^=(t) and ^ = (t)" (5 ' 38) 
By neglecting n 2 = 0(m 2 ), we can approximate (2.44)as 

(ft " XX ] )& = ^ (5-39) 

with 

E l = m 2 or m 2 . (5.40) 
To the same approximation, (2.40) gives 

Ci = -xUi/n- (5.41) 
In terms of the Pauli spin matrices Gi, write 

h - XX t = a + + ^2^2 + (5.42) 
where, neglecting 0(m~ 2 ), 

a = ^ + r?(/ ( 2 + /|) 2 + 2^/2] (5.43) 

and 



6l= T?(f? + f 2 p +%ffif e f p , 

b 2 = ntfefp, (5.44) 

&3= -^ 2 + h 2 [(/e 4 "/ p 4 )-2^]. 



The eigenvalues of (5.42) are 



m„ = a + 6 



and (5.45) 



ml = a — b 
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with 



b = (bl + bt + bj)i 



(5.46) 



In the charged lepton sector, it is convenient to introduce the polar coordi- 
nates ai and fii through 



b\ = b sin ol\ cos Pi, 62 = b sin ai sin and &3 = &cosa/. 

-2> 



(5.47) 



Neglecting 0(m ), we find the components of the eigenfunctions ip e and ip^ 
in ( 5.38) to be 



1 



and 



sin 

-e l @ l COS 



e ^ cos \oti 



(5.48) 



sin 

There are arbitrary phase factors in iV e and iV M , which will be discussed in 
Appendix C. Here we simply set these normalization factors to be 



N e = N^ = l + 0(m- 2 ). 
Combining (5.38) with (5.48), we have 



(5.49) 



< e il e 

< P||e > 

< Pi \e > 



> \ 
I 



sin ^a>i \ 
-e llSl cos \a\ 

C e 



(5.50) 



and 









( e 1/31 cos \oii N 


Vv = 






sin ±07 




V < p l\»> J 







in which 



C e = 



-m 1 x f 



—e 



sin 

*^ cos ^ck/ 



(5.51) 



(5.52) 



c u = - 



m 1 x f 



e cos \oi\ 
sin ^ol\ 
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and 



x - 



Xe 

Xp 



(5.53) 



given by (2.36)-(2.37). Correspondingly, neglecting 0(m 2 ) we find 



( < ejr > ^ 

< Pi\r > 
<P\\r> 



I m 1 Xe \ 
mT x Xp 

1 



(5.54) 



Define the transformation matrix W to be 



W = 



' <e|e| > <^|e| > < r\t { > N 

< e\Pl > < V>\Pl > < r\pi > 
{ < e|P| > < ^|P| > < r|P| > j 



(5.55) 



with ej_, p[, P[ being the same eigenvectors e, p, P of (5.29). Thus, 



< e * 



= w 



' e x 
\ T 1 



(5.56) 



( e ) 




( e \ 




= w ] 


V 


V r J 







and the lepton mapping matrix is given by 



(5.57) 



with (Ui) given by the HPS form (5.1). 

Combining these matrices, we derive the lepton mapping matrix V/_ map 
given in Table 1. To compare with the experimental neutrino mapping matrix, 
there is still a phase convention which will be discussed in Appendix C. 



5.4 Jarlskog Invariant 
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The Jarlskog invariant J\ can be calculated by using (3.29) and replacing 
Uckm by V/_m a p of Table 1. The result is 

* = 6K^y [1 + GO""*)] (5-58) 

with 

= ^/ e /, + (r,/m)/ e / P [// 2 - rf{f t + / 2 )(/ e - / p )(/ e + 2f p )] (5.59) 
In the limit m — ► oo, J/ becomes 

J, = [6(m 2 - m 2 )]- 1 ™^/^ + 0(m- x )]. (5.60) 
In the limit 77 — > 0, we have 

^ = ^(m^/e/p + m;Ve/p) + 0{rf). (5.61) 

The same result can also be derived by using the perturbative expression (3.34) 
and replacing V in (3.35)-(3.40) by V/_ map . 

Remarks 

As mentioned before, the leptonic timeon model predicts two relations be- 
tween the seven observables (5.25)-(5.27). To see how a test can be made, 
we may take the five parameters in the theory to be 

m, fj,, 77, f £ and f p (5.62) 

of (5.23)-(5.24), or the equivalent set 

m, a, b, ai and fii (5.63) 

with a, b given by (5.43), (5.46) and 07, by (5.47). Using (2.50) and setting 
rriL = m T , the mass of the heavy lepton r, we have for the first parameter in 
(5.63) 

m 2 = m 2 + o(7f). (5.64) 

Likewise, from (5.45) 

a=o( m M + m e) ( 5 - 65 ) 
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and 

b= l -(ml-ml). (5.66) 

The two remaining angular parameters a\ and (3i in (5.63), or the equivalent 
parameters f e and f p can then be determined by using 

Ji and any \V t] \] (5.67) 

i.e., the absolute value of any element of the experimentally measured neutrino 
mapping matrix. All other remaining elements of the neutrino mapping matrix 
may serve as part of the test of the timeon model. Further discussions are 
given in Appendix D. 
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Table 1. Lepton Mapping Matrix (neglecting 0(m 2 )), with the parameters 
a = ai and (3 = (3\ given by (5.47); c e , by (5.52), and Xet Xp by (5.53) and 
(2.36)-(2.37). In the limit m — > oo, c e , c M , m^Xe and wT l Xy a " become 0. 
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Appendix A Two Forms of Mass Matrix 
A.l General Formulation 

As in (2.45), the mass matrix M. and its related Hamiltonian H of a Dirac 
field operator \I/ with n-generation components can be written as 

U = ¥M^ {A.l) 

in which 

M = M\ (A2) 

denoting a hermitian matrix. Decompose \I/ into a sum of left-handed and 
right-handed parts: 

^ = C + 1Z (A3) 

with 

£ = l(l + 75 )ilr and ft = I(i_ 75 )^. (A4) 
Correspondingly, (A.l) becomes 

H = tfMiJl + V)M^L. (A5) 

Assume n > 3 and .M to have an imaginary part so that H\sT,C and CP 
violating. 

A different form of an n-generation T and CP violating mass Hamiltonian 
can be written in the form similar to (1.3), also with a Dirac operator ip of n 
components: 

H = ip\<3^A + iFlAlb)ip, (A6) 
where Q and T are both n-dimensional hermitian matrices, 

g = g^ and T = T\ {A.l) 

For n > 3, Q and T both nonzero, the Hamiltonian H is T, P and CP 
violating. As in (A.3)-(A.4), we resolve ip in a similar form: 

i) = L + R (A8) 
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with 

L = i(l + 75 )^ and fl = i(l- 75 W>. (A9) 
Thus, (A. 6) becomes 

H = L\Q - iT) lA R + R\Q + %T)^L, (A10) 

different from (A. 5). 

In the standard model, excluding the mass Hamiltonian, only the left hand 
components of \ and \ quarks are linked by their VK-interaction. Hence, the 
right-hand component 1Z or R can undergo an independent arbitrary unitary 
transformation. Because of this freedom, we can bring (A. 10) into the form 
(A. 5), or vice versa, as is well known. We shall review this equivalence, and 
then discuss how this equivalence can be altered by imposing new restrictions 
on these matrices. 

To show this, we begin with the form (A. 10). Define 

M = Q-iT (All) 

and assume it to be nonsingular (i.e., the eigenvalues of M^M are all nonzero.) 
On account of (A. 7), the hermitian conjugate of M is 

M* = g + iF. (A.12) 

Since MM^ is hermitian, there exists a unitary matrix Vl that can diagonalize 
MM\ with 

= rn 2 D = Diagonal. (A13) 

For every eigenvector of MM^ with eigenvalue A, the corresponding vector 
M^cj) is an eigenvector of M^M with the same eigenvalue A. Thus, M^M can 
also be diagonalized by another unitary matrix Vr as 

V ] R M^MV R = m 2 D , (A14) 

with m 2 D the same diagonal matrix of (A. 13). 

Multiply (A. 13) on the right by m^ 1 , it follows that 

VlMV R = m D , (A15) 
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provided that we define 

V R = M^V L m- D \ (A16) 

One can readily see that Vl and Vr thus defined satisfies VjVl = 1- VrVr — 1 
as well as (A. 13) and (A. 14). Since R can be transformed independently from 
L, we can transform the ip field by 

L V L L {All) 

and 

R -> Vr#. (A18) 

Next let us examine the mass matrix M. of (A.1)-(A.2). Because M. is 
hermitian, it can be diagonalized by a single unitary transformation V, with the 
left-handed and right-handed components of the field operator \I/ undergoing 
the same transformation; i.e., in contrast to (A.17)-(A.18), we have 

C VT, (A19) 

?e -> vn (A.20) 

and correspondingly 

W -> ^^74^ (A21) 

with being the corresponding diagonal matrix. So far as the mass matrices 
are concerned, we regard these two mass Hamiltonians 7i and H as equivalent, 
if the diagonal matrix mp of (A. 21) has the same set of eigenvalues as those 
in (A. 15). In this case, we can without loss of generality set 

M 2 = (Q - iT){G + %T) (A22) 

and 

V = V L ; (A.23) 

hence (A. 13) becomes 

V ] M 2 V = m 2 D (A24) 

and therefore 

= m D . (A25) 
(Note that M ^ M or M\ even though M 2 = MMK) 
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A. 2 Restricted Class with Q and T Real 

We will now discuss theories in which both matrices Q and T are real; i.e., 

Q = g* and T = T\ (A26) 

Since Q and T are also hermitian; they must both be symmetric matrices. 
In n-dimension, each of these matrices can carry ^n(n + 1) independent real 
parameters, giving a total of n(n + 1) real parameters. On the other hand, 
A4 being a single hermitian matrix consists of only n 2 real parameters. Thus, 
knowing Q and T , by using (A. 22), we can always determine uniquely the cor- 
responding A4, but not the converse, by expanding in power series as follows. 
Decompose the hermitian A4 into its real and imaginary parts: 

M = R + iI. (A27) 

with R and / both real; hence, R is symmetric and / antisymmetric. On 
account of (A. 26), the real part of (A. 22) is 

R 2 - I 2 = g 2 + T 2 , (A28) 

and the imaginary part is 

{R, I} = [G, T\. (A29) 

In what follows, we assume that Q and T are both known, as in the case when 
the mass matrix is given by (1.3). In addition, T can be regarded as small 
compared to Q. Hence, we can expend R and I in powers of T. Write 

R = G + R 2 + R± + Re + • ■ • (A30) 

and 

/ = h + h + h + • • • , (A31) 

with R n and I m to be of the order of T n and T m respectively. Eqs.(A.28) and 
(A.29) give 

{g, h} = [g, n 
{g, R2} = f 2 + il 
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{g, h} = -{R2, hh 04-32) 

{G, R4} = {h, h} - Rl etc. 

As noted before, so far as these mass matrices are concerned, the two for- 
malisms (A.l) and (A. 6) are regarded as equivalent to each other, provided 
that (A. 28) and (A. 29) hold. Then (A. 32) gives the conditions determining 
the series expansions (A.30)-(A.31) of R and / in terms of g and T . 

It is convenient to write g in terms of its eigenvalues v, /i, m and their 
corresponding eigenvectors e, p, P (as in (2.14)): 

g = vee + fipp + mPP. (A. 33) 

Let A be a vector whose k th component is given by the (i, j) th component of 
the commutator between g and T\ 

= djkAk {AM) 

with Cijk = ±1 depending on (ijk) being an even or odd permutation of 
(1,2,3), and otherwise. From (A. 33) and (2.18)-(2.19), we can readily 
verify that 

A k = r{u(e ■ f)(e x /) + M (p • f)(p x f) + m(P ■ f)(P x f))^ (AM) 

Since / is antisymmetric, so is I\. Write its (ij)th component as 

(h)ij = e ijk J k . (A.S6) 

By using the first equation of (A. 32) with (A.33)-(A.36), it can be readily 
verified that J is related to A by 

J • e = (fi + m)~ A • e 

J-p=(m + vY l A-p (A37) 

and 

J ■ P = {u + !i)- l A-P. 
In the same way, we can solve for R2, h, • • •. 
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Appendix B Proof of Eq.(3.34) for Jarlskog Invariant 
B.l Definitions, Corollaries and Conventions 



Let V = (Via) be a 3 x 3 real orthogonal matrix with positive determinant, 
and indices i and a = 1, 2 and 3; hence, 



V = V\ V- X = V and \V\ = l. 
Given any value of %, define i', i" by 

cm" = 1, 



(B.l) 



(B.2) 



so that (i, i', i") is a cyclic permutation of (1, 2, 3). We shall use the same 
definition for each of such similar indices j, k, I, a, (3, 7. Thus, for a given 
j, the corresponding j' and j" satisfy 



ejj>j» - 1, 



;s.3) 



and likewise 

tkk'k" = €U'l" = e aa 'a" = €0(3' 0" = 6 7 y 7 " = 1. 

Furthermore, for any pair i, a, we have by the expression for V 1 



Vi' a > Vi'a" 
Vi" a > Vi"a" 



= \v\(y- 1 ) ai = v ia 



BA) 



on account of (B.l). This identity will enable us to reduce certain quartic 
products of Vi a to triple products of Vi a , as we shall see. 

B.2 Jarlskog Invariant 

Similar to the relation between V of (3.32) and 

U = Uckm (5-5) 
of (3.33), we define W = (W ia ) through 

U = V + iT q W {B.6) 
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where W has the form 



W ia = £ FijV ja - £ FapV ifl . (B.7) 

For our purpose here, it is necessary to specify only that F and T are real and 
symmetric; i.e., 

Fij = Fji = F*j (B.8) 

and 

We also define for any particular pair of indices k and 7, 

J = U kj U k ,yU* ky U^ (5.10) 

and 

Jo = V kl V k ^V k \,V k %. (5.11) 
Thus, substituting (B.6) into (B.10), we find, to first order in r q , 

J-J = ir q A kl (5.12) 

where 

A fc7 = J (— + ^ Wk7 ' Wk ' 1 )- (5-13) 

V Vfc 7 Vfc/y Vfcy Vfc' 7 / 

Note that /c, 7 are subject to the cyclic convention typified by (B.2). [It will 
turn out that A kl is independent of the choice of k and 7, even though this 
is not assumed.] 

By substituting (B.7) into (B.13), we must obtain an expression of the form 

A fc7 = £ A v ,F lv + £ Ayi'Txx (5.14) 

l X 

where the A's and A's are to be determined. From (B.7), (B.13) and (B.14), 
we see that each A is made of terms having the form JoVj a /Vi a - Consider 
A k : in (B.14) we must put I" = k; hence /, I' are k' ', /c" in some order. Thus 
in (B.7), i is either k! or But the index k" does not occur in (B.13). 
Therefore, we have i = k' , j = k" and a = 7 or 7'. Therefore, 



A k = Jo O + Tf^-O 



Vfc" 7 ' T4" 



'7 
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on account of (B.4). Likewise, 

Aw = Vk'jVk'yVk'y 



(5.15) 



(5.16) 



For the calculation is different, even though the result will be similar. 
Note that in (B.14) /, /' can be k, k! in either order, so that in (B.7) i and 
j can also be k, k! in either order. Thus, we now have four terms instead of 
two: 

'Vfc' 7 , Vfcy Vfc'y Vk- 



Aw — Jo 

= Vk> 7 Vk>y 



+ 



(5.17) 



Vfc 7 Vfc/y Vfcy Vfcy 

T4/ 7 T4/ 7 ' ,yy i Vky Vkj 

Vfc 7 Vfcy 7 7 Vfc/y Vfc/ 7 

= —Vk'jVk>yVk"y> — Vfc 7 VfcyT4 7 " = +Vk» 1 Vk»yVk»y 

since V is orthogonal. 

Similar formulas are obtained for A 7 , A' 7 , A" with a change of sign. Since 
(B.15)-(B.17) all have the same form, the choice of k, 7 in (B.10) and (B.ll) 
is immaterial and we have 

J = J + ir q A 

where 



A = E^'Il^a-E^AA'n^ 



iX" 



(5.18) 
(5.19) 



a 



B.3 Applications to quarks 

By consulting (3.7)-(3.15) and (3.17)-(3.25), we find that 



F — 



fifj sj~ faf(3 



rrii + m. 



Ha + M/3 



(5.20) 



where 



with 



/1 = (/e)t> h = (/p)tj /3 = (fp)h 

h = (Ze)i' /2 = ifp)h h = (/p)i 

mi = m u = 0, ra 2 = m c , ra 3 = m t , 



(5.21) 
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Hi = m d = 0, /i 2 = m s , 112, = m b . (5.22) 
Then (3.34) is seen to be (B.14) if we identify 

A C = A 3 , A t = A 2 , Bi = A u 

A s = A 3 , A b = A 2 and B i = A x . (5.23) 

With these identifications, (3.38)-(3.40) are (B.15)-(B.17) and (3.35)-(3.37) 
are the corresponding formulas for A\, A 2 and ^.3. 



Appendix C Phase Convention in V\. 



map 



In order to compare the lepton mapping matrix V/_ map given in Table 1 and 
the experimentally measured neutrino mapping matrix 

U v = (Uij), (C.l) 

there are certain phase conventions. In definition of 

1 < v\\e > < v 2 \e > < z^3 1 e > ^ 
Vi- map = < villi > <v 2 \li> <v-i\n> (C.2) 
v < V\\t > < v 2 \t > < v$\t > 1 

given by Table 1, we made arbitrary phase choices of these statevectors. This 
allows us to consider the following transformations: 



\e> \ 
\/i> 
\t> j 



\e> \ 
\/i > 
\r> j 



(C.3) 



and 



where 



( h > N 


\ 




\ 


W2 > 




h > 


(CA) 


{ h > , 


1 


I h > 


1 


i 




\ 











(C.5) 


\ 
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and 

/ e^ 1 \ 

Q, v = e^ 2 . (C.6) 

v e^ 3 , 

The experimentally measured neutrino mapping matrix U v is related to V/_ map 
by 

*/„ = JW-mapft; 1 , (C.7) 

with the conditions 

\U V \ = 1 (C.8) 
and the following four matrix elements of U y , 



Uu, U12, U23 and C/ 33 



(C.9) 



all real and positive. Thus, (C.8) and (C.9) determine five of the six phase 
factors in (C.5)-(C6). The remaining one 



& + L + £r + m + m + m 



(C.10) 



does not appear in L^. 



Appendix D Test for Leptonic System 



Following the discussion given in Remarks at the end of Section 5, we may 
consider 

Ji and \V n \ (D.\) 
as an example of (5.67), and two other members, say 



|Vi 2 | and |y 2 3 1 



From Table 1 and noting that c e and c M are 0(m ), we have 

. l2 2 2 a 1 2 a 2 a a 
Vii = - sin — | — cos — I — sin — cos — cos p 
11 3 26 2 3 2 2 M 



+ \Re 



c*(2sin- + e l/3 cos 



(D.2) 
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|T . l2 1 . ,a 1 oQ; 2. a a 
V12 = t: sin — + - cos - -- sin- cos — cos p 
11 3 23 232 2 



-\ R e 



c* (sin cos-) (£3) 



and 



| V23 1 =-sm — + Re[c^sm-\. 



2 2 M 2 

in which a = ai, (3 = j3i and all 0(m~ 2 ) terms are not included. Using 
(5.52)-(5.53) and noting that, to our approximation, (2.36)-(2.37) yield 

therefore, in accordance with (5.52)-(5.53) 

c : = ^/p(/ e sinf-/,e-^cosf) 
and (D.b) 

cJ = ^/p(/ e c </J cosf + / P sinf) 

where 

r = 77. (D.6) 
From (D.3) and (D.5) and eliminating /p by (2.21), we find 

5 1 1 

IV11I 2 = cos a H — sinojcos/3 

12 4 3 

- ~(1 - /e " /p)*(/e + 2/ P ) sin a sin 
X 1 

|Fi 2 | 2 = - - - sin a cos (3 

lT c2 ,2 ' 



- 7T-a - # " #) s (/e " /p) sinasin/3 (D.7) 

o 77 / 



and 



IT/ |2 1 1 

IV23I =4-4 COSa 



~— (1 - 7 e 2 - ffflfe a sin f3, 
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in which a = ai and (3 = Pi are defined by (5.47). Thus, |Vii| 2 , | V12 1 2 and 
IV23I 2 are all functions of m, /i, r, f e and f p . 

Finally, the Jarlskog invariant is given by (5.58)-(5.59). These expressions 
provide the explicit forms for the four observables in (D.1)-(D.2). Together 
with (5.45) and (5.64) for the masses of e, fi and r, we have seven observables 
in terms of five parameters. 




Figure 1. A schematic drawing of the quark mass matrix My± = Gyp^ + 
2^7475. The vibration of t(x) is timeon 
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X 



c 




Figure 2. Geometric Representation of Harrison, Perkins, Scott Transforma- 
tion (5.1) and (5.6)-(5.8) 



The axes OX, OY and OZ = -OD represent e, \i and r. An a = 45° left- 
hand rotation (Vj.)o along e takes r to and fi to OA. A second /? = sin -1 
left-hand rotation (V|)o along z/ 3 takes e to z/i and 0^4 to z/2 which is along 
OB. The hidden bases are l\ = e, ^ || 0^4 and = ^3. 
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